fig2a_data = readtable('fig2a.csv');
fig2b_data = readtable('fig2b.csv');

figure; set(gcf,'position',[360,250,560,360]);

subplot(2,1,1);
plot(fig2a_data{:,1},fig2a_data{:,2},...
    'o','MarkerSize',4.5,'linew',1,'color',[0.000 0.447 0.741],'MarkerFaceColor',[0.0745,0.6235,1]);
set(gca,'fontsize',13);
xlim(600*[-1,1]);
grid on;
set(gca,'xticklabel',[],'ytick',0:0.5:2);
hold all;
plot(fig2a_data{:,3},fig2a_data{:,4},...
    '-','linew',2,'color',[126,47,142]/255);
plot(fig2a_data{:,5},fig2a_data{:,6},':k','linew',3);
legend({'experiment','simplified model','numerical model'},'position',[0.6399    0.7640    0.2958    0.1806]);
ylim([0,2.25]);
ylabel('$|1-\mathcal{C}|$','Interpreter','latex');
set(gca,'position',[0.1268,0.5838,0.8274,0.34]);

subplot(2,1,2);
plot(fig2b_data{:,1},fig2b_data{:,2},...
    'o','MarkerSize',4.5,'linew',1,'MarkerFaceColor',[0.0745,0.6235,1]);
set(gca,'fontsize',13);
xlim(600*[-1,1]);
grid on;
hold all;
plot(fig2b_data{:,3},fig2b_data{:,4},...
    '-','linew',2,'color',[126,47,142]/255);
plot(fig2b_data{:,5},fig2b_data{:,6},':k','linew',3)
plot(269.4357*[1,1],2*pi*[-1,1],':k','linew',3);
ylabel('$\mathrm{arg}(1-\mathcal{C})$','Interpreter','latex');
xlabel('\omega, two-photon detuning [Hz]');
set(gca,'ytick',pi*(-1:0.25:1),'yticklabel',{'-\pi','-3\pi/4','-\pi/2','-\pi/4','0','\pi/4','\pi/2','3\pi/4','\pi'});
set(gca,'position',[0.1268,0.1485,0.8274,0.34]);

ylim([-2,1]);
text(-1400*600/1100,5.1+0,'\bfa','fontsize',14);
text(-1400*600/1100,1.25+0,'\bfb','fontsize',14);
